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Abstract. Photoionisation is a process where absorption of one or several photons 
liberates an electron and creates a hole in a quantum system, such as an atom or a 
molecule. Is it faster to remove an electron using one or many photons, and how to 
define this time? Here we introduce a clock that allows us to define ionisation time for 
both one-photon and many-photon ionisation regimes. The clock uses the interaction 
of the electron or hole spin with the magnetic field created by their orbital motion, 
known as the spin-orbit interaction. The angle of spin precession in the magnetic field 
records time. We use the combination of analytical theory and ab-initio calculations 
to show how ionisation delay depends on the number of absorbed photons, how it 
appears in the experiment and what electron dynamics it signifies. In particular, we 
apply our method to calculate the derived time delays in tunnelling regime of strong- 
field ionisation. 
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1. Introduction 


Quantum mechanics uses operators to predict an outcome of a measurement. However, 
some figures of merit, e.g. the phase of quantum electromagnetic held or the times 
of quantum transitions are not associated with operators. Their measurement has to 
rely on the operational approach, i.e. a measurement protocol yielding a particular 
observable at the detector. Measuring ionisation times with attosecond accuracy has 
been the focus of several recent experiments [l]-6|, using different measurement protocols 
in different ionisation regimes. 

For one-photon ionisation the definition of ionisation time Tws, referred to as 
Eisenbud-Wigner-Smith time, is established 
of recent experiments 


and verified in the analysis 15,9f 12 


. This definition links T\y$ to the phase of the photoelectron 
wave-function </> through its derivative with respect to electron energy: Tws = —df/dE. 

Establishing such a link between the classical concept of time and the parameters of a 
quantum wave-function in the regime of strong field ionisation is the goal of this paper. 

To define the ionisation time for one-photon and multiphoton ionisation regimes 
within the same protocol, we extend the idea of the Larmor clock, originally introduced 
in [l3] to define the time it takes an electron to tunnel through a barrier (see e.g. 14-L6]). 


The Larmor clock measures rotation of the electron spin in an external homogeneous 
magnetic field acting exclusively inside the barrier. The angle of rotation is the hand 
of the clock. Here, we introduce the analogue of this clock, which is based on the 
spin-orbit interaction and is naturally built into many atoms. Physically, the spin-orbit 
interaction can be understood by considering an electron with angular momentum l 
orbiting around the nucleus. In the reference frame associated with the electron, the 
nucleus rotates around it and creates a current. The current creates the magnetic field. 
Precession of the electron spin in this field records time. 

2. Spin-orbit Larmor clock for one-photon ionisation: calibration of the 
clock 


To illustrate how the clock works, consider a Gedanken one-photon ionisation 
experiment, where an s-electron is removed from an atom (e.g. Cs) by a right circularly 
polarised light field, see Fig. [l](a) . There is no spin-orbit interaction in the initial state- 
the ground state of a Cs atom. Spin-orbit interaction turns on upon photon absorption, 
since the electron angular momentum changes from l = 0 (s-state) to l = 1 with its 
projection on the laser propagation direction mi = 1. Thus, photon absorption turns 
on the Larmor clock: the electron spin starts to precess. Spin-orbit interaction is short- 
range, localised near the core. Once the freed electron leaves this area, spin precession 
stops and the clock turns off. 


However, the original Larmor clock 13 used homogeneous magnetic held and hence 


the clock hand rotated with constant speed. In our case, the spin-orbit interaction is 
inhomogeneous, requiring calibration of the clock: the mapping between the angle of 
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Figure 1. Gedanken experiment for calibrating the spin-orbit Larmor clock 
in one-photon ionisation. (a) Cartoon of the experiment, for a Cs atom: a 
circularly polarised field removes an electron from the valence 5s shell, prepared in 
a coherent superposition of the spin-up and spin-down states. In the final state, 
the electron spin has rotated by angle A <j>so- (b) The spin-orbit clock operating 
on the electron viewed as an interferometer, a simplified view. The two arms 
correspond to the spin-up and spin-down pathways, with phase difference A 4>SO] 
Formally, the interferometer describes the following interfering pathways, where 
mj = 1, m s = m' s = ±1/2 correspond to two different arms of the interferometer: 

R 


'i,rn' s ,f REmi = 0,m s ,g) cx ( m' v m' sl f 


mi = 1, m s ,g). Here R is the radial 


part of the dipole operator, S is the angular part of the dipole operator, g , / are 
the radial parts of the initial and final state wavefunction. We have used that 
S|?7p = 0,m s ) oc |to; = l,m s ). (c) Detailed view of the spin-orbit interferometer. 
The spin-down path is itself a double arm, since the spin-down electron (m s = —1/2) 
and the final orbital momentum To; = 1 can proceed via both j = 1/2 and j = 3/2 
continua. The single (spin-up) arm and the double (spin-down) arm interfere in 
the final continuum state with To; = 1. Formally, it corresponds to the following 


interfering pathways: R, = 1 ,m s ,g), where m[ = 1 

m s = m! s = ±1/2, correspond to two different arms of the interferometer. 


rotation and the ionisation time. To calibrate the clock we consider rotation of electron 
spin during one-photon ionisation, where the ionisation time Tws is well established. 

Let us prepare the electron in the initial spin-polarised state, |si n ) = (1 /a/2)[| — 
1 / 2 ) + 11 / 2 )], with the phase (p characterising the initial orientation of the spin 

in the polarisation plane, and calculate the angle of rotation of the electron spin in 
this plane during ionisation. The final state, for the final orbital momentum mi = 1, 
is |sfi n ) = (1 /\/ 2 )[ 0 , 4 ,| — 1 / 2 ) + 11 / 2 )], where the and are the ionisation 

amplitudes for the spin-up and spin-down ionisation pathways. The spin has rotated by 
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the angle A (ft so = argfa^a^], equal to the phase delay between the spin-up and spin-down 
ionisation pathways. 

The amplitudes cq, cq, have been originally derived by U. Fano [l7j: cq = R 3 , 
cq = |(i ? 3 + 2 i?i), where i?i i3 are the radial transition matrix elements into the 
degenerate continuum states with total angular momentum j = 1/2 and j = 3/2. 

To find the phase difference we need to find arg[cqcr(]: 


a t a l 


Rsl 3 (^ + 2 ^) 


= 3 (W 2 + 2 |i?i||i? 3 |(cos( 0 f - <j>?) + isin(</>f - c/>f))), ( 1 ) 


A fso = arg R 3 
Equation (( 2 ]) yields: 


3(^ + 2^ 


tanA0 5 o = 


= arctan 


2 i?i 

w 

sin^^ — ) 

\Rs\ 

2 + 2 

|i? 3 |(cos( 0 f - 0 f) 


sin(—A 0 i 3 ) 


( 2 ) 


( 3 ) 


0.5|i? 3 |/|i?i| + cos(0i 3 ) 

where the phase difference A</>i 3 is defined as A0i 3 = (f)f — the relative phase 
between i? lj 3 (E). Their dependence on electron energy E is very similar, up to a 
small off-set A Eso due to the spin-orbit interaction in the ionisation channel: R 3 (E) = 
Ri(E — AE S o ) [l7|. The phases are shifted accordingly: f 3 (E) = <j)f(E — A E S o) §] 
Using Taylor expansion, we find A</i 3 = (E) — </>f \E — A Eso) — ~ T wsAEso, where 
Tws = —d(j>i/dE is the Wigner-Smith ionisation time 17,8,11 . Thus, we have connected 


the angle of spin rotation during ionisation A fso to the Wigner-Smith ionisation time, 
calibrating our clock, 

sin^wsAUso) 


tan A f>so = 


( 4 ) 


0.51 7? 3 1 /1 | + cos(t ws AE so ) ' 

The inhomogeneous character of the spin-orbit interaction makes the relationship 
between A <fso and Tws nonlinear, introducing the extra term 0.5|i? 3 |/|i?i| in the 
denominator, but does not invalidate the clock. 

How does the spin-orbit interaction measure the ionisation time? The spin-orbit 
clock works like an interferometer, see Fig. 0 b, c). The angle of rotation is given by 
the relative phase between the spin-up (parallel to the orbital momentum) and spin- 
down pathways. The spin-up pathway proceeds only via the j = 3/2 continuum. The 
spin-down pathway is a double arm: it can proceed via both the j = 3/2 and j = 1/2 
states. The spin-orbit interaction introduces the phase delay A0 33 in the double arm. 
This small perturbation records the ionisation time, with Eq. Q connecting A (f>so an d 
TWS- 

f The photoionisation matrix elements we discuss here are complex, and we deal with their phases. 
The phase lag in a real-valued radial wave-function, corresponding to a given ionisation channel j = 3/2 
or j = 1/2 discussed in 17 , translates into the phase of the complex-valued photoionisation matrix 


element, leading to the phase difference A0 13 of the corresponding photoionisation matrix elements 
discussed here. 
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Figure 2. Gedanken experiment for measuring ionisation time in strong-field 
ionisation, (a) Cartoon of an experiment, for a Kr atom: Intense circularly polarised 
field removes a p~ electron from the valence 4 p 6 shell via absorption of many photons, 
creating a rotating hole, (b) The spin-orbit clock operating on the hole as an 
interferometer, for an electron removed with mi = —1 and the hole left with Ml = — 1. 
As in Fig. EK the spin-up path is a double arm, since for Ml = —1 the spin-up 
pathway (Ms = 1/2) can proceed via both J = 1/2 and J = 3/2 hole states. The 
ionisation amplitudes, up to the angular coefficients relating the orbital momentum 
L , spin S, and the total angular momentum J, are (for J = 3/2) and T-f" (for 
J = 1/2). The relevant angular (Clebsch-Gordan) coefficients for each pathway are 
also indicated separately. 


3. Strong-field ionisation 


We now turn to strong-field ionisation in intense IR fields-the regime of recent 
experiments |3,61 aimed at measuring ionisation times and focus on the definition of 
the ionisation time in this regime. 

To address this problem, consider ionisation of a Kr atom [Fig. |2][a)] . The ground 
state of Kr + is spin-orbit split by the energy A Eso = 0.665eV into the P 3 / 2 and Pi/ 2 
states with total angular momentum 3/2,1/2. Both will be coherently populated by 


ionisation 18 in a few-cycle intense IR laser pulse, with coherence approaching 90% for 


nearly single-cycle pulses 18 . The loss of coherence arises when photo-electron spectra, 
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Figure 3. Strong-field ionisation of Krypton by a a single-cycle, circularly polarised IR 
pulse. Angle and momentum-resolved photo-electron spectra calculated numerically 
for the two ionisation potentials corresponding to the two ground states of Kr + , -P3/2 
(I p = 14.0 eV) and P 1 / 2 (I p = 14.67 eV). Red solid contours correspond to P3/2, 
blue dashed contours correspond to T\ / 2 • The inner-most contour corresponds to 
0.9 level, other contours shown in steps of 0.1. The pulse had a vector-potential 
A lR) = —(F/w) cos 4 (wt/4)(coswf x + sinwty) with F = 0.05 a.u. and w = 0.057 a.u. 
The radial coordinate gives the electron momentum in atomic units. 


correlated to two different core states, do not overlap or overlap only partially. 

Figure [3] shows nearly complete overlap of these spectra for the ultrashort circularly 
polarised pulse used in our ab-initio calculations, see Appendix A for the numerical 
details, confirming nearly 100% coherence of the hole motion in this case. Note that the 
lack of 100 % coherence affects the amount of the coherently moving charge in the ion, 
but not the timing of its dynamics. Thus, even for coherences below 100%, one is still 
able to use the rotation of the spin of the hole, triggered by the spin-orbit dynamics of 
the electron charge in the ion, as a clock. 

At the same time, the spin-orbit interaction in the ionisation channels becomes 
completely negligible in strong fields. The importance of this effect can be gauged using 


the Analytical A-matrix Approach (AAM) (T9 23 , which has been verified against ab- 
initio simulations in 


23[24|. Application of ARM to calculating ionisation phases is 


described in |Appendix B In the tunnelling limit, AAM yields the following expression 


for the relative phase between two degenerate continuum states associated with electron 

2.3 x 10” 7 rad, 


total momentum j — l + 1/2 and j' = l — 1/2: £so ~ 0.2la 2 F 2 /1^ 2 


r\j 
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which is completely negligible. Here a is the fine structure constant. For the estimates 
we used typical values of ionisation potential I p ~ 0.5 a.u. and the strength of the laser 
field F = 0.06 a.u. 


see 


Appendix B for the details of the derivation). 


The formal description of the spin-orbit interferometer in Fig. [2] is similar to the 
one-photon case in Fig. [TJ There is no spin-orbit interaction in the ground state of 
Kr: the filled valence 4p 6 -shcll has equal number of p~ and p + electrons ‘rotating’ 
in opposite directions. Ionisation by a nearly single-cycle, right-circularly polarised 
IR pulse breaks the balance between the co-rotating and counter-rotating electrons: 


intense right-circularly polarised IR pulse preferentially removes the p electron 25 -27 
and creates a p - hole. This starts the clock. The angle of rotation of the hole spin at a 


time-delay r after the IR pulse is (see Appendix C): 


tanA0so = 


sin (A E S oT ~ A0i 3 ) 


0.5|T 3 I/ITj I + cos(A E SO T - A 01 3 ) 


( 5 ) 


The dependence of the ionisation dynamics on the IR pulse intensity, duration, shape, 
etc, is fully encoded in the matrix elements Tf and Tf [Fig. |2](b)] describing strong- 
field ionisation amplitudes for the removal of mi = —1 electron, leaving the hole 
in pj,Mj valence spin-orbitals with total angular momentum J = 3/2 and J = 1/2 
correspondingly, for a given final electron momentum p at the detector. It corresponds 
to population of P 3/2 and P 1/2 states of Kr + . The phases of T 3 R are </> 3 and 
A 0 i 3 = <j>J — (f) 3 . For the p + electron (mi = 1) the expression is similar, except that Tf 3 
are different. The tiny phase shift between spin-down and spin-up ionisation amplitudes 
correlated to the P 3 /2 state of Kr + due to spin-orbit interaction in ionisation channels 
has the same estimate as above, £so 0.21 a 2 F 2 /Ip /2 ~ 2.3 x 10 'rad, and is also 

negligible. 

The clock angles in Eqs. (|3]), ([5]) are virtually identical, except for the term A E so t 


describing the hole dynamics 28 upon ionisation. The analogy in angle-time mapping 


in Eqs. (J3]) and (J5]) allows us to establish the definition of strong-field ionisation time. 
Indeed, Eq. ([3]) calibrates the clock and establishes the mapping between the angle of 
spin rotation and the ionisation time. Eq. ([5]) contains the same mapping. Thus, the 
time of hole formation is encoded in A 0 13 = (f>J—f 3 , accumulated in the second (double) 
arm of the interferometer. We shall now analyse these phases to extract the strong-filed 
ionisation time. 

The phases <fr[, (jff encode the electron interaction with the potentials of the 
core states Pl / 2 and P 3 / 2 . These potentials have two contributions, = U c +Uf 3 . Here 
U c is common for both states and is dominated by the long-range Coulomb potential, 
while Uf 3 are different for the two core states, reflecting different spatial distributions 
of their electron densities, see Appendix D Thus, fj 3 = £> 33 + (ftf 3 . 


For the same final kinetic momentum p of the continuum electron, in the strong- 
field ionisation regime, the difference between and <f 3 comes from slightly different 
ionisation potentials into the Pi/ 2 and P 3 / 2 states: (j) 3 = f c (I p ) and <f 3 = f> c (I p — AE S o), 
see Appendix B Hence A <f\ 3 = — 0 3 ~ A Esod<f c /dl p and one should convert A0 33 
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into time, dividing by the energy of spin-orbit splitting [note the ‘minus’ sign in Eq. (|5])]: 

= A( P c i3 = d( P c 

TSI AE so dip 

Equation (J6]) defines ionisation time in the strong-field regime. 

The second part of the relative phase, Aff 3 = — results from the different core 

potentials for the P 3 / 2 and Pi / 2 core states, e.g. due to the different angular structure 
of the electron density. It does not depend on A E S ch he. the period of the clock, and 
hence can not be converted into the time-delay in the formation of the hole. 

We note that the derivation presented here is not applicable in the weak-field regime. 
Firstly, the neglected spin-orbit interaction in ionisation channel may become important 
in weak fields. However, it is a rather standard approximation to ignore spin-orbit 
interaction in the ionisation channel compared to the spin-orbit interaction for the core 
electrons (see e.g. |29|), since the core electrons are plenty and stay near the core, where 
the spin-orbit interaction is strong, while the sole continuum electron leaves the core 
region. Secondly, the explicit dependence of phases </>i j3 on ionisation potential, used 
in deriving Eq. ([ 6 ]), arises naturally only in the strong-field regime, see 
Therefore, the expression Eq. (J 6 ]) may not hold in the weak field regime. 


Appendix B 


4. Reading spin-orbit Larmor clock in strong-field regime 

In the Gedanken experiment described in Section [3j Eq. (J5| is sufficient to introduce 
the strong-field ionisation time Tsi by comparing it to Eq. Q which has calibrated 
the clock. In contrast to the one-photon case, where the clock stops as soon as the 
liberated electron leaves the range of the spin-orbit interaction potential, the discussion 
in Section [3] does not involve stopping the clock. Indeed, the clock operates on the 
hole states J = 3/2 and J = 1/2. The hole spin periodically rotates after ionisation is 
completed. Thus, the clock continues to count time after it has recorded the rotation 
angle related to the ionisation time. To stop the clock and read the information out we 
can apply the second pulse. It allows us to get direct access to the phase A0 13 , which 
records time in the spin-orbit interferometer shown in Fig. [2](b) . 

Here is how it works. Consider the pump-probe experimental scheme shown in 
Fig. 0 a,b). The pump, which starts the clock, is a nearly single cycle right circularly 
polarised IR pulse. The probe, which stops the clock, is a left circularly polarised 
attosecond XUV pulse. It comes with an attosecond-controllcd delay r and fills the 
hole in p~ orbital by exciting an electron from a deeper s-orbital (see Fig.4(a)). More 
generally, the probe pulse promotes the core into an excited S'-state, where the spin- 
orbit splitting is absent, e.g. 4s4p 6 or any other suitable state. Broad bandwidth 
of the attosecond probe pulse couples both P 3 / 2 and Pi / 2 to the same final S'-state, 
as in [18] (see Fig. 4(b)). As opposed to the Gedanken experiment above, in 
laboratory experiments the initial spin-up and spin-down components of the ground 
state are incoherent, and the single arm of the interferometer yields background for the 
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Figure 4. Schematic of a laboratory pump-probe experiment implementing the 
spin-orbit Larmor clock for strong-held ionisation, for Kr atom, (a) Cartoon of the 
experiment. Multiphoton ionisation with a strong, right-circularly polarised infrared 
pump pulse creates a p“-hole and starts the clock. Attosecond extreme ultraviolet 
probe fills the p _ -hole by promoting an electron from the inner s-shell. This transition 
stops the spin-orbit clock, since spin-orbit interaction for s-states is absent, (b) 
Analysis of the experiment as a two-path interferometer. Two pathways via the P \/2 
and P 3/2 states of the ion interfere in the final S-state of the Kr + . 


interference in the double arm in Fig. [f](b). Left-polarised probe ensures that the final 
S'-state can be reached only if the electron missing in the Kr core after ionisation is the 
p~ electron. 

The population w of the final state is (see Appendix E): 


w ~ l^i| 2 + 1^4.3 1 2 + 2 1^4i11^4.31 cos(AE i s'07‘ — A0i3(p)) + |A31 


(7) 


Here, A 1 = 2T, (p)d 1/2 F u> (Qi)^2/27 and A 3 = T 3 (p )d 3 / 2 F u (Q 3 )y/2j27 are the 
transition amplitudes for the two interfering pathways corresponding to the removal of 
the spin-down p~ electron. I11 addition to the multi-photon ionisation matrix elements 
Tf 3 (p) they include the real-valued radial transition matrix elements d \/2 and d 3 / 2 
between the P1/2, T3/2 and the final S'-state of Kr + , and the spectral amplitudes of 
the attosecond pulse, F u (...), at the excitation energies Q 3) i from the P 3 /2,i/2 states to 
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the final 5-state. The background |A 3 | 2 , A 3 = T 3 _ (p)d 3 / 2 F w (fi 3 ) 1/2/27, corresponds to 
the removal of the spin-up p~ electron (See Appendix E for details of the derivation). 
Modulation of w versus the pump-probe delay r yields the phase A0 i3 (p) between A\ 
and A 3 . It can be measured, e.g. using attosecond transient absorption (l8|. The 
sensitivity of the phase A0 13 (p) to the final momentum p, and the errors that it can 
introduce into the transient absorption measurement of the ionization time, are discussed 
below. 

The phase A0i 3 includes two contributions: (i) the relative phase due to the same 
core potential in both ionisation pathways, A0 33 , which can be translated into time- 
delay and (ii) the phase A0 33 , related to the different electron-core potentials in the 
two ionisation channels. This phase reflects correlation between the electron and the 
core and can not be translated into time. If an experiment does not distinguish between 
these two contributions to A0i 3 , the phase A0f 3 related to the electron-hole correlation 
will look like a time shift. 

Figure |5|a) shows how the total phase (red squares), which can be measured by 
transient absorption, and its two separate parts A</> 33 (blue circles) and A</>f 3 (green 
triangles), depend on the laser wavelength, i.e. the minimum number of photons 
N = Ip/oj required to reach the ionisation threshold, for fixed laser intensity. Figure[5](b) 
shows Tsi (blue circles) and apparent time delays r eh = — A0f 3 /A E S o (green triangles). 
The apparent delay T e h is not negligible for P 3 / 2 and Pi/ 2 states of Kr + . 

To obtain results in Fig. [5](a,b), we have calculated the phases accumulated due 
to the Coulomb potential and the short range components of the core potential for the 
two ionisation channels, corresponding to the ionic states P 3 / 2 and Pi/ 2 . Note that the 
short-range potentials in these two channels are different, see Appendix D To obtain 
time-delays, we have divided the relative phases by the difference in the ionisation 
potentials, A E S o- The phases were calculated using the APM method [l9]-|23], for the 
characteristic momentum of the photo-electron distribution po = Aqsivl)i(ujtt)/(ojtt), 
where A 0 is the amplitude of the field vector potential and tt = Im[p(p 0 )] is the so-called 


‘Keldysh tunnelling time’, the imaginary part of the saddle point t s (p 0 ), see Appendix 


[C| For this momentum, which is very close to the peak of the distribution for the short 
laser pulse, the ionisation phases have simple analytical expressions in the tunnelling 
limit: 


A« 3 ^ -A Eso/r/, (8) 

A</4 ~ —OAF 2 /Ip/ 2 . (9) 

Note that A0^ 3 is proportional to A Eso and therefore leads to proper time-delay, while 
the phase A0 33 accumulated due to the different short-range potentials does not scale 
with A Eso and cannot be translated into proper time. Thus, every time is phase, but 
not every phase is time. 

Since transient absorption experiments do not detect the final energy (or 
the momentum) of the electron, we have also checked that the phases and the 
resulting times are only very weakly sensitive to the final electron momenta within 
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Figure 5. Analysis of time delays in strong-field ionisation, (a) Calculated phase 
A(f>i 3 for the pump-probe experiment, as a function of the minimum number of photons 
required for ionisation, N = I p /ui, I p is ionisation potential, ui is laser frequency. The 
calculations were done for a Kr atom and the circular field intensity 2.5 x 10 14 W/cm 2 . 
Blue circles show phase associated with the actual time-delay. Green triangles show the 
phase that does not correspond to time-delays but is a leftover from the electron-hole 
correlation. Total phase is shown as red squares, (b) Real (blue circles) and ’apparent’ 
(green triangles) ionisation delays as a function of the number of photons required 
for ionisation, N. (c,d) Physical picture underlying the results: iV-dependence of the 
electron exit position from the potential well (see Appendix B) (c) and the cartoon of 
the ionisation process (d). 


the region surrounding the peak of the photo-electron signal. This analysis is 
presented in Fig. [6j where the ionisation time-delays are overlayed with the electron 
spectrum generated by the single-cycle pump pulse with the vector-potential A L (t) = 
— A 0 cos 4 (cuf/4)(cos(u;t) x + sin(u;f)y), with A 0 = F/u , F = 0.05 a.u. and u> = 0.0465 
a.u. 

The difference in ionisation times within the full width at the half-maximum of 
the distribution is ±5 asec. This number provides an estimate for possible errors in 
transient absorption measurements of ionisation delays caused by averaging over the 
photo-electron distribution. Note that such measurements will also inevitably include 
the apparent delays r e h = —Af>f 3 /AEso associated with the phase Aff 3 . For the 
specific example shown in Fig. [5])b) r e h ~ 10 asec. Importantly, in the tunnelling limit 
Tsi is intensity-independent while the apparent delay r e h is proportional to the laser 
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Figure 6. Dependence of ionisation delays tsi on the electron momentum 
at the detector, for Hydrogen atom. The inner-most contour in the electron 
spectrum corresponds to 0.9 level, other contours are shown in steps of 0.1. The 
color bar shows tsi in attoseconds. The pulse had a vector-potential A^(f) = 
— (F/w) cos 4 (wf/4)(cos(wf) x + sin(wt) y) with F = 0.05 a.u. and w = 0.0465 a.u. 
Results are obtained using the ARM theory. 


intensity, see Eq. (J9]) . This factor might be used to separate these two contributions. 

Results presented in Fig. |5](b) show that, as we increase the laser wavelength A and 
hence the number of photons N = I p /u required for ionisation, the ionisation time in 
Fig. [5|b) decreases. This dependence has simple explanation. As A decreases, the laser 
frequency uj increases, ionisation becomes less adiabatic and the electron splashes out of 
the potential well closer to the core, see Fig. [5j^c,d) . From there, it runs to the detector, 
accumulating the phase and consequently the time-delay tsi■ The closer the electron is 
launched, the larger is the accumulated phase. Note that no delay is accumulated under 
the barrier, see Section [6] for details. 
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5. Attoclock measurements of strong-field ionisation delay 


The spin-orbit Larmor clock has offered us a general procedure for defining ionisation 
times in both one-photon and strong field ionisation regimes. Using the same general 
procedure we have arrived at two different expressions in the weak-field one-photon 
ionisation regime and in the strong-field regime. In the weak field regime we found the 
Wigner-Smith ionisation time. In the strong-field regime we found an expression which 
agrees with the result of a completely different derivation described in 23 . Importantly, 


while we have derived ionisation times using spin-orbit interaction, our results do not 
depend on it. Therefore, the detection of the strong-field ionisation time does not have 
to rely on the spin-orbit interaction. 

Consider, for example, the so-called attoclock setup 13,30 , which measures angle- 
resolved electron spectra produced in nearly circular, few-femtosecond IR pulses. Such 
pulses send electrons released at different instants of time in different directions, 
providing the link between the direction of electron velocity at the detector and the time 
of its release. Nearly single-cycle pulse creates preferred direction of electron escape, 
from which the ionisation delay can be reconstructed [3,23,30 


The angle 


at 


which the majority of electrons are detected, relative to the detection angle expected 
in the absence of the core potential, is called the off-set angle. We now show that </> max 
can measure the ionisation delay tsi derived above, provided that effects leading to 


transient population trapping of released electrons and ’’negative” ionisation times 23 


are negligible. In particular, such regime can be achieved in the long wave-length limit 
(but is not limited to it). 

To this end, we consider the benchmark system-the hydrogen atom, where fully 
ab-initio simulation of ionisation dynamics in the strong circularly polarised IR field is 
possible. We solve the time-dependent Scrodinger equation numerically exactly and use 
results of the numerical experiment to find 0 max . Details of the calculation are described 
Appendix A Red circles in Fig. [7] show the ionisation delay At = (</> max — 59)/< 


m 


L0 


extracted from the ab-initio photoelectron spectra, where u is the laser frequency and 
the small correction 59 to the off-set angle f) max is introduced by the rapidly changing 
pulse envelope of the nearly single-cycle laser pulse we have used 23 . The blue curve 


with squares, which shows Tsi Eq. i§. lies on top of the ab-initio results. To calculate 
tsi analytically, we have used the ARM theory |23j. 

We stress that the definition of Tsi is not restricted to the long-wavelengths limit 
shown in Fig. [7} It is only the ability of the attoclock set-up to measure exclusively 
this time delay, without additional contributions associated with transient population 
trapping in Rydberg states leading to negative ionisation times 23 , that has restricted 
our consideration to the wavelength regime shown in Fig. [T} Nevertheless, it is important 
to demonstrate at least one example, where the time delay tsi derived from the idea of 
the spin-orbit Larmor clock can be experimentally or numerically detected. 
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Figure 7. Attoclock measurements of strong-field ionisation delay, (a) Red circles 
show numerically calculated envelope-free offset angle 4>max — 56. The small correction 
59 to the off-set angle 4>max is introduced by the rapidly changing pulse envelope of the 
nearly single-cycle laser pulse 23 and is subtracted from (j>max to present envelope- 
free results for the offset angle. The blue squares connected by the blue line show 
u>\tsi\, the offset angle corresponding to time-delay |rs/|, to is the laser frequency, 
(b) Red circles show the ionisation delay At = ( (f>max — 59) /uj extracted from the 
ab-initio photoelectron spectra. The blue squares connected by the blue line show 
r si ■ All calculations were done for a hydrogen atom and the circular field intensity 
1.75 x 10 14 W/cm 2 . 


6. Strong-field ionisation delay and tunnelling delay 


Strong-field ionisation is often viewed as tunnelling through the barrier created by the 
binding potential and the laser electric field. While our analysis has never relied on 
the tunnelling picture, our definition is consistent with the Larmor time —d<p/dV for 
tunnelling through a static barrier of height V 14,15 , equal to I p in our case (see 
Fig. §a ,b)), f> is the phase of electron wave-function. However, Fig. [8|a,b) emphasises 
the difference in the two processes, which is in the boundary or initial conditions for the 
tunnelling dynamics. In Fig. [8](a) , the electron current is incident on the barrier and it 
can lead to the appearance of tunnelling delays, i.e. phase and time delays accumulated 
during the motion under the barrier. In Fig. |8](b) the tunnelling starts from the real- 
valued wave-function of the bound state. It is a plausible assumption that in this case 
and for the low-frequency laser field, the polarized bound state carries negligible current 
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Figure 8. Cartoon illustrating the analogy and the difference between (a) the standard 
barrier penetration problem, and (b) optical tunnelling through the barrier created by 
the laser field and the core potential in strong-field ionisation. 


incident on the barrier, and that tunnelling would occur from the tail of the initial wave- 
function already present under the barrier. Since the initial wave-function is real-valued 
in the barrier region, the phase in Eq. ([3]) may get no contribution from the tunnelling 
region, leading to no delay associated with the under-barrier part of the electron motion. 
Indeed, the analytical calculation of the phase in Eq. ([3]) yields no contribution from 
the under-barrier region, at least in the regime of Fig. [7} 

As it follows from excellent agreement between analytical and numerical results in 
Fig. 0a ,b), the analytical calculation of the phase is accurate, and optical tunnelling is 
not associated with time delay. The delay Tsi is only due to electron interaction with 
long-range core (Coulomb) potential and is explicitly accumulated after the exit from 
the barrier. 


7. Conclusions 


We have illustrated the concept and the meaning of time delays in strong held ionisation. 
In one electron systems, these delays are related to electron interaction with the nucleus. 
In the tunnelling limit, comparison of numerical and analytical results unambiguously 
demonstrates the absence of tunnelling delays. Non-equilibrium charge dynamics 
excited in a many electron atom or a molecule by the laser held and the electron- 


electron correlations 20,31 could lead to additional phase 5f> 32 and additional delays 


Srsi = —ddf/dlp contributing to Tsi■ Our work shows why and how ionisation delays 
provide a window into such dynamics in complex systems. 

Production of a coherent superposition of many ionic states and hence of coherent 
hole dynamics is the key aspect of interaction with ultra-short light pulses. Any pump- 
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probe experiment resolving these dynamics aims to find phases between the coherently 
populated states. As a result of electron-core correlations, not all phases are mapped 
into time: the formation of the hole wavepacket is characterised not only by the overall 
time-delays, but also by additional phases accumulated during the ionisation process 
due to the different core potentials for the different final states of the ion. 

What do these phases mean? Given that the electron wavepackets correlated to 
different core states overlap at the detector, the hole presents a coherent wavepacket 
characterised by the relative phases of its different spectral components. Analysis of 
spectral phase is common in characterisation of ultrashort pulses in optics. Linear 
spectral phase records the arrival time, while non-linear phase is associated with pulse 
dispersion. Such dispersion is the closest analogue of the phase shifts related to electron- 
core correlations. 
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Appendix A. Ab-initio calculations 


The numerical procedure and the code are described in detail in 23,33 . The method 


has been monitored for convergence by changing the maximum angular momentum up 
to L m ax = 300, the radial grid size was increased up to r max = 2500 a.u., and by varying 
the step size of the radial grid Sr down to 0.05 a.u. In the presented calculations, the 
step size of the radial grid was Sr = 0.15 a.u., the time-step was St = 0.04 a.u., the box 
size was 1500 a.u., and L max = 150. 

For Hydrogen atom, the spectrum was obtained by projection on the exact field- 
free continuum states of the //-atom after the end of the laser pulse. The photoelectron 
spectra include the volume element oc jf, both in numerical and analytical calculations. 
The volume element shifts the position of the peak of the distribution and thus 
affects the off-set angle, however, in the exact same way for numerical and analytical 
spectra. In these both numerical and analytical calculations we define the laser field 
Fl(£) = —<9A L{t)/dt via the vector-potential A ^(t): 


A L (t) = — A 0 f(t) (cos (uit) x + sin (art) y), 
where f(t) is the pulse envelope and uj is the carrier frequency, 


(A.l) 


f(t ) = cos 4 (a;t/4). 


(A.2) 
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For Kr atom, the calculations have been performed using the effective one-electron 
model potential 


U Kr (r) = 


1 + (36 — 1) * exp(— rjr) 


+ u 0 , 


(A.3) 

We follow the recipe 
35 , using the additional tuning potential U 0 which is added only at 


based on the DFT potential used by D. Bauer and co-workers 34 
described in Ref. 


the first radial grid point ri = 0.5 (the radial grid step was Ar = 0.05 a.u.) and is equal 
to zero everywhere else. The parameter r/ = 2.64343586965 a.u. has been adjusted 
to yield the correct ionisation potential of Kr for the lowest J = 3/2 ionic state, with 
additional fine-tuning achieved by setting Uq = 0.0249a.u., giving I p = —0.5145022731 
a.u. For the J = 1/2 core state the tuning potential was adjusted to U 0 = 22.7629 a.u., 
yielding I p = —0.5389895221 a.u. 

The photoelectron spectrum was calculated by propagating the wavefunction for 
sufficiently long time after the end of the laser pulse (typically 2 cycles, the convergence 
has been monitored up to 10 cycles), then applying a spatial mask to filter out the central 
part of the wavepacket within 100 Bohr from the origin, and performing the Fourier 
transform of the remaining part of the wavepacket. We have independently validated this 
procedure using the Hydrogen atom, where it has been calibrated against the projection 
of the wavefunction on the exact scattering continuum states for Hydrogen. The mask 
radius was chosen based on this calibration in Hydrogen. 

Appendix B. Calculation of the phase accumulated due to interactions in 
ionisation channels 

Appendix B.l. Definition of the strong-field ionisation phase accumulated due to 
interactions in ionisation channels 

To evaluate the relative phase between the two ionisation channels in Kr, we use the 


R-matrix based method (ARM) 119,20 generalised for the case of circularly polarised 


fields 21-23 


The ARM method allows one to obtain an analytical expression for the total phase 
accumulated in each ionisation channel: 


f>j(P,t s (P,I P )) = 


' t s —iK 


dtUji / dCv(C) 


(B.l) 


where Uj( r) is the potential defining the interaction, n = \J2I P , I p is the ionisation 
potential, v(t) = p + A(t), A(f) is vector-potential of the laser field, T —* oo is the 
observation time and p is the electron final momentum at the observation time. 

The time t s (p,/ p ) (see 


231) is the complex-valued solution of the saddle point 


21 


equations for the ionisation in circularly polarised field: 


dSy(T, p, t s ) 
dt x 


= /, 


pi 


(B.2) 
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where Sy (T. p) is the Volkov phase accumulated by the electron in the laser held only: 


S v (T,p,t s ) = ^ f dt[p + A{t))' 2 . 

J to 


(B.3) 


The coordinate of exit presented in Fig. |5](c) of the main text is: 

/•Re[t s ] 

d(v( (). 


ro = 


(B.4) 


Since t s (p,I p ) depends on I p , the phase (j)j(p, t s (p, I p )) also depends on I p . The phase 
difference in the two channels is accumulated due to the different If s: the difference in 
ionisation potentials leads to slightly different t s and thus slightly different trajectories in 


the two channels. These trajectories are the arguments of Uj in Eq.(B.l). The common 


part of the phase is accumulated due to the Coulomb potential. The channel-specific 
part is accumulated due to the channel-specific core potential discussed in Appendix 
|Dl The phase accumulated due to spiu-orbit interaction in the ionisation channel is 


negligible and is estimated below in Appendix B.2 


Appendix B.2. The phase accumulated due to spin-orbit interaction in ionisation 
channel 


We estimate the relative phase between the two ionisation channels corresponding to 
spin-up and spin-down ionisation pathways, with orbital momentum l and two values 
of electron total momentum: j = 1 + 1/2, and j' = l — 1/2. We use the spin-orbit 
interaction potential: 


V S o(r ) = 


+ -1(1 + 1) -s(s + l) 

4c 2 r 3 


(B.5) 


where c = 1/a ~ 137, in atomic units, s = 1/2 is electron spin (a being the fine- 
structure constant). The phase difference fso is expressed via the difference between 
the potentials corresponding to j and j': 

l + 1/2 


A V so (r) = 


2 c 2 r 3 


(B.6) 


We now calculate the phase difference using Eq. (B.l) connecting the phase to the 
potential. Substituting the electron trajectory in the tunnelling limit r = r 0 + Ft 2 j 2, 
where r^ = I p /F, F is the held strength, we obtain the following integral: 


(l + 1/2) r°° dt 
2c 2 r 3 J 0 (1 + g) 3 ’ 


(B.7) 


where l, the electron angular momentum along the trajectory, remains constant in the 
pure tunnelling limit. Evaluating the integral: 



dt 


(i + gf) 


r dx 

F J„ (1 - x‘y ' F ’ 


(B-8) 
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we obtain: 


fso — 


(l + 1/2) J2T p (l + 1/2) F 2 

V ' -0.59^—^ = -0.42 v J ’ 


2 c 2 rg 


5/2 ' 


(B.9) 


Note that in general the angular momentum of the electron l is changing with time and 
should be included in the integrand. However, the integral is accumulated in the vicinity 
of the core and therefore for estimates in the tunnelling limit we can use the value of 
angular momentum Iq when the electron exits the tunnelling barrier. In the tunnelling 
limit Iq —>• 0, since when the electron exits the tunnelling barrier its velocity is parallel 
to electron displacement. Thus, for typical field strength F = 0.06 a.u. and I p = 0.5 
a.u. the phase difference £so ~ 0.21F 2 /(c 2 I p ^ 2 ) ~ 2.3 x 10 -7 rad. and is completely 
negligible. 


Appendix C. Rotation of the hole spin in strong field ionisation: Gedanken 
experiment in Kr atom 

Consider Kr atom in its ground state. There is no spin-orbit interaction in the ground 
state of the neutral Kr: the P-shell is filled by 6 p-electrons, with equal number of 
p~ and p + electrons ‘rotating’ in opposite directions. Ionisation by strong, circularly 
polarised IR laser field breaks the balance between p~ and p + electrons 25 and starts 
the spin-orbit Larmor clock. Intense right-circularly polarised IR pulse prefers to 
remove the p~ electron |25j 26], i.e. mi = —1. Let us set the initial spin state 
to be |s in ) = ot\ — 1/2) + e^/3 11/2), where a,/3 are real numbers and the phase </> 
characterises the initial orientation of the spin. Once the p~ electron is removed, 
the quantum state of the core acquires uncompensated angular momentum, with the 
hole created with Ml = — 1 and uncompensated spin. The spin state of the hole is 
|s in ) = a\- 1/2) + e^d|l/2), since the spins and the angular momenta of the hole and 
the electron are the same at the moment of separation. As this state is not an eigenstate 
of the Hamiltonian, the hole spin starts to precess. 

We shall now calculate the angle of rotation of the hole spin. The final spin 
state for the fixed orientation of the final orbital momentum Ml = — 1 is |sfi n ) = 
a±a | — 1/2) + <qe^/3|l/2), where the a^ and cq are the strong field ionisation amplitudes 
for the spin-up and spin-down ionisation pathways. We first specify our notations and 
introduce the ionisation amplitude T~(I P ) corresponding to the removal of p~ electron in 
the absence of the spin-orbit splitting of the core state. Here I p is the ionisation potential. 
The amplitudes T 3 and T 1; which include the spin-orbit splitting, are proportional to 
T~(I p ): Tf oc T~(I p ) and Tfi oc T~(I p + A E S o), and they correspond to the removal 
of the p~ electron 125,28 . 

Full ionisation amplitudes into the hole states Pj,Mj include the projections 
(LM l , SM s \ JMj) given by the Clebsch-Gordan coefficients, C JI f Ij ,,, , with M L = 1. 
Taking these projections into account, we find that the amplitude of ionisation into the 
hole state J = 3/2, Mj = 3/2 is Tf ■ The amplitude of ionisation into the state J = 3/2, 
Mj = 1/2 is the amplitude of ionisation into the state J = 1/2, Mj = 1/2 is 
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— y |T 3 _ . Now, we project these states back onto the L, Ml, S, Ms basis to find cq 
and cq. This yields the amplitude to find the hole angular momentum Ml = — 1 and 
Mg = —1/2 at time t cq = Tfe~ lE3 ^ 2t , while the amplitude to find the hole angular 
momentum M L = —1 and M s = 1/2 at the time t is op- = | (2Tf e~ lEl / 2t + Tf e~ lE3 / 2t ). 
Here £3/2 is the energy of the ground state, \Ei/ 2 \ = \E 3 / 2 \ + AEso- 

To establish the rotation angle (the phase of the double arm of the interferometer 
relative to the phase of the single arm (see Fig. 2(b)); the single arm now corresponds 
to spin-down pathway, while the double arm now corresponds to spin-up pathway) we 
need to find argfaqal]: 

a t a t = 3 (1^3 | 2 + 21^1 || r 3 |(cos(03 —4*1 + 

AE so t) +i sin(</»3 - 4>i + AE so t)), (C.l) 


arg 


t 3 i jpr +2 t , 




= arctan 


2 \ T i ll r 3 I sm(<j% - + AE so t) 

|T 3 -|2 + 2|T 1 ~||T 3 “|(cos(03 - + AE so t) 


Equation (C.2) yields: 

tan A 4>so = 


sin(A E so t - Af 13 ) 


0.5|T 3 |/|T\ I + cos(A£ 5 o^ - A0i 3 ) ’ 
where the phase difference A0 i3 is defined as A0i 3 = — f> E . 


(C.2) 


(C.3) 


Appendix D. Core potentials in two different ionisation channels 


To illustrate the effect of electron-hole correlations on definition and measurement of 
time, we consider the contribution of the channel specific core potential Vljmj { r ), that 
arises from the Coulomb interaction between the electron and the core. This potential 
has the following form: 


V LJ Mj( r) = / dr 


Pti 


r — r 


= / dr 


' || r — r r || ( eJ Mj\ r ')( r '\ eJ Mj) = 


ftJMj / 2l 71 /[' 

’lm' l ,\m' s 'lm l ,\m s \2 s 

T' A/T„ AJ' ' 


m l ,m l Ms,m s 

OO 


CXJ nj 1 

J2 p Li( cos/3) / 
Li=0 Ldo 


dr r 


I 12 


L\ 
pL l+l 


-M s 
2 6 


\R(eJMj-r')\ 2 + I drr 


dSlYl M ,(e',4f)YL Mh {& A ')x 


/ ,/2 


r* L ^ 


/\ |2 


pf L\ -|-1 


\R(eJMj-r')\ 


(D.l) 


where, L — J ± 1/2, is the orbital angular momentum hxed for a given spin-orbital, ft 
is the solid angle between the vectors r and r', and can be written as cos (3 — r • r', e 
represents the effective principle quantum number corresponding to the energy of the 
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spin-orbital under consideration, and R(eJMj\r ) is the radial part of the wavefunction 
associated to the said spin-orbital. 

Including all terms together, we have: 


v LJ mao = y czumC j lZum s y 


An 


m s 


--- 2Ll + 1 

Li= 0 M Ll =—Li 


E YI iMli (6^)(R Li )x 


dSlY’ m (#,&). (D.2) 


Here <f is the angle in polarisation plane, 6 is the angle calculated from the laser 
propagation direction, and (Rh) is the expectation value of the radial component, as 
calculated using the Roothaan-Hartree-Fock (RHF) orbitals, defined as: 


r l { t ) = ^2 c i p c iq 

p,q 


1 

ry*L~\-\ 


7((«i, + K P r ) + rLr (( K v + K i,) r ) • 


(D.3) 


where, c ip , c iq are the coefficients for the Slater-Type Orbitals (STO) and i p , i q the 
corresponding indices defining the nodes in the wavefunction under consideration, used 
for the RHF calculations [36 , and 7 is the lower, whereas T is the upper incomplete- 
gamma function. Taking into account Wigner 3j-coefficients from the integral: 


dnYi MI (e'A')Y LML {e'^')Y LlMLi { cose') = (-i) M ^(2L + i) 


2L\ + 1 


47T 


■X 


L Ll L\(L L x L\ 
M l M Li -M' l ) \ 0 0 0)' 


(D.4) 


we obtain the selection rules. For Li, the selection rules are: (a) 2 L + L x is even (so 
only L 1 even are allowed in the summation over L\) and (b) the triangle inequality 
\L — L\\ < L < L + L\ which gives 0 < L\ < 2L. For all other cases the integral is zero, 
and M' l = M l + M Ll ■ 

Taking L\ = 2L', the expression for Vljmj is: 


Vi 


LJMj 


r = 


L 2 L' 

( 2 ^+ 1 ) E E E 

Ml,L '=0 —— 2L' 

L 2II L 


fiJMj 

LM l +M 2L ,±MsLM l ,\M s 


An 


AL' + 1 l M l M 2 L' 


-M' 


L 21! L 
0 0 0 


YwM, r AW)(Rw)- (D.5) 


From the Clebsch-Gordan coefficients, we have two conditions on Ml and Ms for a 
given Mj\ 


Ml + M 2 L' + Mg — Mj, 
Ml + M s = Mj, 


(D. 6 ) 

(D.7) 
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which can only be possible if Mil' = 0. The final expression is: 


V L jMj(r) — (2 L + 1) 


L 

£ £ 

Ml,M g L '—0 


(-i) 


m l 


(jJMj 


LM l , r m s 


47T 


L 2 L' 


4L' + 1 \ M l 0 -M l 
L 21/ L 


x 


0 0 0 


Using the definition of Y 2 l'q, we can simplify further to give: 


\Y* m (9A)(R 2L ,) (D.8) 


Vljmj (r) — (2 L + 1) £ £(-d 

M L ,M s L ’=0 


M l 


qJMj 


LM l A M s 


2 I L 2 L' 


M l 0 -M l 
L 21/ L 


x 


0 0 0 


P 2L '(cos 9) (R 2 l')- (D.9) 


Note that L' = 0 corresponds to Coulomb potential, common in both channels. Consider 
the case when ionisation liberates the p + electron (L = 1) populating the hole states 
J = 3/2,1/2 and Mj = 1/2 (the result for p~ is the same). For the calculation of the 
difference between two core potentials we use the same trajectory with averaged I p . The 
corrections associated with the difference in the trajectories are of higher order and are 
not included here. 

The difference in core potentials for this trajectory is: 


^ 1 , 3 / 2 , 1/2 ( r ) ~ ^ l , l / 2 , l / 2 ( r ) — 3 ^ 


L '=0 


1 2 L' 1 
0 0 0 


P 2 L'(cos6)(R 2L ')x 


.Ml,M s 


^ 3/2 1/2 

1 m l ,\m s 



£ (-i) M ‘ 

M l , M s 


^ 1 / 21/2 

i 1 

21! 

1 )] 

1 M L ,\M S 

\M l 

0 

-Ml) 


(D.10) 


As expected for the common Coulomb potential, the difference for II = 0 is zero: 


£ (-D"‘ 

M L ,M s 


( 1 ° 

1 n 

^3/2 1/2 

2 

^1/21/2 

2 " 

l M l 0 

-M l / [ 

1 M L ,\M S 


1 M L ,\M S 




+ (-l) 


1 0 

1 0 -1 



0. (D.ll) 
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The only term left is the one corresponding to II = 1, which gives 


n,3/2,./2(r) - n,iAi/ 2 (r) = 3| J 2 J ]f> 2 (cos0)(fl 2 }x 


Ml, Ms 


L 2 
M l 0 


1 

-M l 


)\ 

^3/2 1/2 

2 

a4/21/2 

2" 


1 m l ,\m s 


1 M l ,\M s 




3a/— P 2 {cos9){R 2 ) 



2 1 
153 + 


-l)\ — 
ni 30 


= le 2 (cos#)(fl 2 ) = -Iff 


since for 6 = 7 t/2 , P 2 = —1/2. The expression for (R 2 ) is: 


12 r °o v -| / -| 

JJ 2' |n/,r.J\|2 , / ill 2' n/,r. |2 1 / /4\r , 2 ' 1 


(D. 12 ) 


« 2 = / AV"-|fl(6L;r')r+ / *V“-|fl(eL;r')|" = + rg - ) , (D.13) 

Jo r Jr r r \ r / r 

which can be found from the incomplete gamma functions. The difference between the 
two core potentials is: Vi 3 / 2 1/2 (r) — Vl 1 / 21/2 (r) ~ —4.444/(10r 3 ), since (R 2 ) = 4.444 a.u. 


for Kr 37 and the contribution of the second term in Eq. (D.13) vanishes for r —> 00 . 
To calculate the respective relative phase A(j)f 3 , we use Eq. (B.l) and substitute the 
difference in short range core potentials given above. 


Appendix E. Pump-probe signal: the details of derivation 


The goal of this section is to derive population in the final E-state of the Kr ion at the 
end of the pump-probe experiment, see Eq. Q of the main text. 

For a laboratory experiment, we need two requirements. First, we want to turn 
on and turn off the clock on demand, i.e. we need to stop the rotation of the core 
spin on demand. Second, we would like to measure the phase A0i 3 directly. The 
second condition is satisfied automatically, since the initial superposition of spin up 
and spin down states is incoherent and therefore the single arm of the interferometer 
(in Fig. [2](b) ) will not interfere with the double arm in a real experiment. Thus, the 
laboratory experiment will only record the interference in the double arm, and the 
single arm will give background. To start the clock, we apply a nearly single-cycle 
right circularly polarised IR pulse to create a p-hole. To stop the clock, we apply a left 
circularly polarised laser field to induce a transition from the s-shcll of the Kr atom, 
filling the Ml = — 1 hole in the p-shcll and leaving the hole in s-orbital. There is no 
angular momentum in the s-hole, and there is no SO interaction. Thus, the left-circular 
probe stops the clock that was started by the right circular pump. 

For a fixed final state of the continuum electron, characterised by momentum p at 
the detector, the population S = \ui/2,m l =o\ + |c r -i/ 2 ,M i =o| in the final s-state can be 
obtained using the following equation: 


°’m s ,m l =0 


* (WanM d 


(E.l) 
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where T; on = J2jm, a JM J 'f’jM J e~ lEjt is the coherent superposition of the two ionic 
states, created after ionisation, for a given final momentum of the electron at the 
detector. Here ajMj is the complex amplitude of ionisation into hole state 
The wavefunction v hfi n (t) = represents the final S'-state of the ion. 

Following Fig. 4(b), we write one particle spin-orbitals participating in the one- 
electron transition from the inner valence s-orbital to the outer-valence J, Mj orbital 
ifjMj and fjs as a product of angular and radial wave-functions: ipjMj = 
ips — '4>s{r)\LMi n SMs ). Taking into account that L — 0, Ml = 0, S = 1/2 in the final 
state, we obtain: ips — fAs( r )|0 0,1/2 Ms). 

The dipole operator can be factorised into the radial and angular parts, d = r 
where £ = 1 corresponds to the right polarised pulse, £ = — 1 corresponds to the left 
polarised pulse: = d x + ifd y . Evaluating the integral over t, we rewrite the equation 

in equivalent form: 


M s ,M l -0 — F u {Es — Ej)ajMjdj 

x (JMj\l M L ,1/2M' s ) 

x (lM L ,l/ 2 Ms|SJ 00 , 1/2 M g 


(E-2) 


Here, F w (Es — Ej) is the Fourier image of the probe pulse F asec (f) taken at the transition 
frequency. The real-valued radial matrix element dj = (ipj(r)\r\ips( r )) describes the 
electron transition from the inner-valence state S to the final sate | JMj), leaving the 
hole in S'-orbital. For the left circularly polarised held, the angular part of the dipole 
operator ^1 M L = — 1,1/2 M' s S_! 0 0,1/2 M$j = Sm s ,m' s VN , N = 2/3. Thus, we 
obtain 


° r M s ,M L = 0 — £ VNdjF u (Es — Ej)ajMj(JMj\l —1,1/2 Ms). (E.3) 

J,Mj 

The Clebsch-Gordan coefficients C/f s=1 / 2Ms = (JMj\l —1,1/2 Ms) are equal 

, r 3/2-3/2 _ , ^3/2-1/2 _ , ,fo ^1/2-1/2 _ _ 

LU ’ L 1,-1,1/2 -1/2 — L 1,-1,1/2,1/2 — 1 /VO, Oi ~ 

Left polarised pulse promotes s-electron to Ml = —1 hole, right polarised pulse 
to Ml = 1 hole, linearly polarised pulse to Ml = 0 hole. Thus, if right circularly 
polarised pulse is used as a pump, left circularly polarised probe will probe ionisation of 
p~ electron, whereas right circularly polarised probe will probe ionisation of p + electron. 
Linearly polarised probe will probe both p + and p~ pathways at the same time. For the 
left-circularly polarised probe we obtain: 



2 

<* 1 / 2,0 


N 

9 ^ 3/2 


—dl 


Tf | F u {Q 3 )f 


(E.4) 
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and 

k-iAoP = ^4 2 |T s -| 2 |F„(Si 3 )| 2 + ^<i? /2 |Tf| 2 |F„(n 1 )| 2 + 

AN . 

—^1/2^3/21^ || T 3 |F w (fi 3 )||F w (fi 1 )| cos(AE so r - A0i 3 ), (E.5) 

where r is the time of arrival of the attosecond pulse, = E s — E 3 / 2 , ^1 = E s — Ey 2 , 
Es is the energy of the final S'-state. Transform limited attosecond pulse is assumed for 
this calculation. The population in S'-state is S — \(J\/2,m l =o\ + |c r -i/2,M I# =o| , and can 
be measured by transient absorption of the XUV probe. 
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